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Abstract 

We briefly outline the models of the nn transition in nuclei and analyse a well known 
calculation based on diagram technique published in 1992 and repeated recently by V. 
Kopeliovich and I. Potashnikova. 
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1. At present several models of nn transitions in medium are treated. The part of them 
gives radically different results. This is because the nn transition in nuclei is extremely sensitive 
to the details of the model and so we focus on the physics of the problem. 

As far back as 1992, the model of nn transitions in light nuclei based on diagram technique 
for direct reactions (later on referred to as the model 1) has been pubhshed [1,2]. In 1996 this 
calculation was repeated for deuteron [3]. However, in [4] we abandoned this model for reasons 
given below. The model based on field-theoretical approach with finite time interval (model 2) 
has been proposed [4-6]. 

"Critical examination of the field-theoretical approach to the neutron-antineutron oscilla- 
tions in nuclei" [7] is that the calculation of [1-3] has been repeated. Also several "objections" 
against [4-6] are adduced. They are wrong or irrelevant (see section 2). For instance, in the 
abstract [7] we read: 

1) "Infrared divergences do not appear within the correct treatment of analytical properties 
of the amplitudes." 

The analytical properties of the amplitudes have nothing to do with the infrared divergences. 
The infrared divergences do not appear since the model 1 used by authors [7] is, by construction, 
infrared- free (sec section 3). The analytical properties of the amplitudes are related to the 
calculation of the integral only (which is trivial). 

In this paper we compare the model based on the field-theoretical approach (section 2) with 
the model based on diagram technique for direct reactions (section 3) and show that the latter 
model is unsuitable for the problem under study. Also it is shown that interpretation of a 
number of important points of the above-mentioned models given in [7] is wrong (sections 2-4). 
Section 5 contains the conclusion. 

2. Model 2 describes both nn transition in the medium [4,5] (Fig. la) and nuclei [8] 
(Fig. lb). It corresponds to the standard formulation of the problem: the \in >-statc is the 
eigenfunctions of unperturbed Hamiltonian. In the case of the process shown in Fig. la, this 
is the neutron plane wave: 

np{x) — exp(— ipx), (1) 

p — (e„,p„), Cn = p^/2m -|- [/„, where Un is the neutron potential. In the case of Fig. lb, 
this is the wave function of bound state [8]. For the nucleus in the initial state we take the 
one-particle shell model. Since the neutron is described by the bound state wave function and 

not the plane wave, the virtual antineutron in Fig. lb is described by Green function and not 
propagator. (The Feynman propagator is the partial case of Green function [8,9].) 
We consider the process shown in Fig. la. The interaction Hamiltonian is 
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(2) 



Here Hnn and "H are the Hamiltonians of the nn transition [10] and the n-medium interaction, 
respectively; e is a small parameter with e = 1/r, where r is the free-space nn oscillation time; 
"ifn and "iffi are the operators of the neutron and antineutron fields; m„ = nin = m. 




Figure 1: nn transition in the medium (a) and nuclei (b) followed by annihilation. 

In the model 2 the antineutron propagator can be bare or dressed. The corresponding results 
differ radically. If the antineutron propagator is bare (the antineutron self-energy E = 0), the 
(S'-matrix amplitudes corresponding to Figs, la and lb contain the infrared singularity. It is 
conditioned by zero momentum transfer in the nn transition vertex. To avoid the infrared 
singularities the problem is formulated on the finite time interval (t, 0) [6]. Then the matrix 
element of evolution operator U{t,0) is calculated (see Eq. (13) of Ref. [4]). The lower limit 
on the free-space nn oscillation time is found to be [4,5] 

t'' = 10'^ yr. (3) 

This value is interpreted as the estimation from above. 

The fact that the evolution operator is used is not new [9,11]. Recall that the most part 
of physical problems is formulated on the finite time interval. The model described above is 
standard and so it created no questions up to now. However, in [7] we read: 

2) "The author of [4-6] tries to reconstruct the space-time picture of the process..." 
In reality we do not more than calculate the matrix element of evolution operator. 

3) "...new rules seem to be proposed in [4,5] instead of well known Feynman rules." 

We calculate the matrix element of evolution operator only. In regard Feynman rules, 
the problem formulation (process model) in the quantum field theory and nuclear physics 
differ principally. In particular, there is density- dependence (coordinate-dependence) in nuclear 
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physics. This is absolutely different problems and so the Feynman rules developed in quantum 
electrodynamics (QED) can be only element of the corresponding model of nuclear reaction or 
decay. This model (more precisely, the model of nn transitions in light nuclei based on diagram 
technique for direct reactions) was published in [1,2] and analyzed in the next section. So the 
words "instead of well known Feynman rules" mislead the reader at least. 

4) " These new rules should allow to reproduce all the well known results of nuclear theory" . 

The consideration of the matrix element of evolution operator is interpreted [7] as "new 
rules". As for the second part of the remark, indeed, all the well known results of nuclear 
theory have been obtained by means of (S-operator ([/-operator). 

The "remarks" 1) - 4) are due to the misunderstanding of the authors of [7]. Also we recall 
that the above-given problem formulation created no questions up to now. If the antineutron 
self-energy equal to zero, the process amplitude is singular. In this case the calculation of matrix 
element is really non-standard [4,5]. However, the pp. l)-4) are irrelevant to the calculation of 
matrix element. 

For the model 2 with dressed propagator the calculation is standard [5,13]. The lower limit 
on the free-space nn oscillation time is 

t'^ = 1.2 • lO'' s. (4) 

The models with bare and dressed propagators were analyzed and compared in [5]. 

3. Model 1 is shown in Fig. 2. Denote: A = (N.Z), B = [N — 1,Z) are the initial and 
intermediate nuclei, M is the amplitude of virtual decay A ^ n+{A — l), M^^'> is the amplitude 
of nB annihilation in (n) mesons. En is the pole neutron binding energy; m, tua, Tns are the 
masses of the nucleon and nuclei A and B, respectively. The process amplitude is [2] 

= / dadE M{q)Mi:^){mA) 

27r4 J ^ (q2 - 2mE - i0)2[q2 + 2mB{E + £;„) - iO)] ' ^ ' 

For deuteron the process probability Wiit) is given by Eqs. (6), (12) and (18) of Ref. [2]: 

^l(i) — ^ d-^mesonst 1 
^ d-^mescms ^ ,-,9 ^npi (6) 

where Vnp is the np annihilation width. This is final result. At the same time the authors [7] 
write "There is no final formula for Td^mesons in [2]." 

Wc list the main drawbacks of the above-given model which are essential for the problem 
under study: 

a) The model does not reproduce the nn transitions in the medium and vacuum. If E^. 0, 
Wi diverges (see also Eqs. (15) and (17) of Ref. 
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Figure 2: Model 1 for the nn transition in the nuclei followed by annihilation. 

b) Contrary to the model 2, the amphtude (3) cannot be obtained from the Hamilto- 
nian because in the interaction Hamiltonian there is no term which induces the virtual decay 
(N,Z) ^n + (N -1,Z). 

c) The model does not contain the infrared singularity for any process including the nn 
transition, whereas it exists for the processes in the medium and vacuum (see [4-6]). This 
brings up the question: Why? The answer is that for the propagator the infrared singularity 
cannot be in principle since the particle is virtual: ^ rr? + p^. Due to this the model is 
infrared-free. On the other, hand the neutron is in the bound state and should be described by 
wave function and not the propagator. 

d) Since the model is formulated in the momentum representation, it does not describe the 
coordinate-dependence, in particular the loss of particle intensity due to absorption. Also there 
is a no the dependence on nuclear density! The model is crude and has very restricted range 
of applicability. 

We consider the points b) and c). The nn transition takes place in the propagator. As the 
result the model is infrared-free. For the processes with zero momentum transfer this fact is 
crucial since it changes the functional structure of the amplitude. 

On the other hand, the neutron propagator arises owing to the vertex of virtual decay 
A — )■ n + (A — 1). However, as pointed out above, in the interaction Hamiltonian there is no 
term which induces the virtual decay A ^ n + {A — 1). This vertex is the artificial element of 
the model. It was introduced in order for the neutron (pole particle) state to be separated. 

Wc assert that for the problem under study this scheme is incorrect. The neutron state is 
described wrongly. The diagram technique for direct reactions has been developed and adapted 
to the direct type reactions. The term "diagram technique for direct reactions" emphasizes this 
circumstance. The processes with non-zero momentum transfer are considerably less sensitive 
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to the description of pole particle state. The approach is very handy, useful and simple since it is 
formulated in the momentum representation. It was applied by us for the calculation of knock- 
out reactions and p-nuclear annihilation [12]. The price of simplicity is that its applicability 
range is restricted. At the same time, as is seen from p. c), the process under study is extremely 
sensitive to the description of neutron state. The same is true for the value of antineutron self 
energy [13]. Besides, the problem is unstable [13]. 

Finally, since the operator T-Lnn acts on the neutron, in the model 1 the vertex of virtual 
decay n + (A — l) is introduced because one should separate out the neutron state. This 
scheme is artificial because in the interaction Hamiltonian there is no term which induces the 
virtual decay A — > n + {A — 1). The neutron of the nucleus is in the bound state and so 
it cannot be described by the propagator. More precisely the model containing the vertex of 
virtual decay is a crude one. It is inapplicable for the problem under study. Alternative method 
is given by the model 2 which does not contain the above-mentioned vertex. The shell model 
used in the \in >-state of the model 2 has no need of a commentary. 

For above-given reasons we abandoned this model [4]. The model based on the potential 
description of n- medium interaction (potential model) has been considered in [5,14] in detail. 

4. The rest of the assertions of [7] should be commented. 

5) The words "we give some general arguments..." are wrong since these arguments result 
from the model 1. 

Figure 1 adduced in [7] determines uniquely the process model (model 1). Equation (8) and 
subsequent consideration relate to the model 1. 

6) The consideration of "energy scale" [7] is nonsense. What is the "energy scale" for nuclear 
/3-decay, for example? 

7) "The reason of suppression is the localization of the neutron inside the nucleus." 

To understand that the localization of the neutron inside the nucleus is unrelated to the 
suppression, it is suffice to compare the nuclear /3-decay and decay of free neutron. The sole 
possible reason discussed in the literature [14,15] is the antineutron potential. In the models 
1 and 2 the potential is not introduced and above-mentioned mechanism of suppression is 
inoperative. We also recall that the nn transition in the nuclei followed by annihilation is the 
dynamical, two-step process with the characteristic time Tch ~ 10~^^ s [4,5,13]. 

8) "If the infrared divergence takes place for the process of nn transitions in nucleus, it 
would take place also for the nucleus form-factor at zero momentum transfer. But it is well 
known not to be the case, as we illustrated in this section" 

So the authors illustrated that the model 1 is infrared-free and the large part of [7] is devoted 
to this "problem" . However, this trivial fact is well known [9,11]: the infrared divergence cannot 
take place for virtual particle. 
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The errors indicated in pp. 1) - 8) and the fact that the main calculation is the special 
case of the model proposed by us [1,2] in 1992 are incompatible with the paper title "Critical 
examination..." and the beginning of the acknowledgements "Present investigation, performed 
partly with pedagogical purposes" . 

5. In conclusion, the process under study is extremely sensitive to the description of neutron 
state (model 1 or 2). The same is true for the value of antineutron self-energy E. If E changes 
in the hmits 10 MeV > E > MeV, the lower hmit on the free-space nn oscillation time r is 
in the range [8,13] 

10^^ yr > T > 1.2 • lO'' s. (7) 

For the free-space ab oscillations the picture is exactly the same: the ab transition probability 
is extremely sensitive to the difference of masses — m^. The fact that process amplitude is 
in the peculiar point (sec Fig. 1) is the basic reason why the range (7) is very wide. The values 
r*^ = 1.2 ■ 10^ s and = 10^® yr are interpreted as the estimations from below (conservative 
limit) and from above, respectively. The estimation from below t'^ = 1.2 ■ 10^ s exceeds the 
restriction given by the Grenoble reactor experiment [16] by a factor of 14 and the lower limit 
given by potential model [17] by a factor of 5. At the same time the range of uncertainty of r 
is too wide. Further theoretical and experimental investigations are desirable. 
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